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ABSTRACT 


«rody„amlc forces 

Se"r“-^ ^~-s’as one o“The1et^“fHr f 
dynamic stabiUty analyses! etltfness coefficients required for rotor- 

high ™ed’’‘‘axLraorte^?" experimental measurements made on a small, 

tn^nvLtJLtf tie'iauJ^tJ ^rMr-tSLa^Sl 


INTRODUCTION 


There are many mechanisms which cause rotordynamic instahlHM^« ^r, ^ >-k- 
and compressors An ^nQf-aK^ i -u-rr u ^*^uyn<uuic xnscapiiities in turbines 

they limit the operation of rnl'e ^"=“‘';J“«s can become a costly problem since 
They can be destructive If the In^Xty^SLllf 

machinery! “D^yU^rrrlsrtrtJ^^Ss'tLre^rLtLirir^h^^^^ 
d1lnrL°‘°'o::f S"tr"°"’ production and heivy m^SnTn1e“\^^^^^^ 

pr^uc't^^f '"rprod^ctL''^'""^ '=»* hy-produ«“^om"c"d; otl 

Problem :uh this‘^c»pres:or i-tablllty 
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In mathematical terms an instability occurs when an eigenvalue of the equation 
that describes the motion of the rotor has a positive real component. In engi- 
neering terms, the nonsynchronous amplitude of rotor whirl (any amplitude which has 
a frequency not equal to running speed) grows with time until the system exceeds 
its operating limits. In many instances the rotor orbits can be large enough to de- 
stroy bearings and seals. 

Computer programs are now available [3] which can predict whirling frequencies 
and threshold speeds, provided that accurate values for the destabilizing coeffi- 
cients are used as input. The destabilizing coefficients usually take the form of 
a linearized cross— coupled stiffness term [4], but sometimes are expressed as neg- 
ative damping [5], depending on the mechanism producing the instability. 

One mechanism, which is a type of self-excited instability, is the aerodynamic 
excitation due to circumferential variation of blade clearance in axial-flow com- 
pressors and axial-flow turbines. An objective of the work reported here was to 
quantify a term in the mathematical equation due to Alford which predicts the mag- 
nitude of the aerodynamic cross-coupled stiffness, along with verifying the form of 
the equation. 

Alford actually presented two mechanisms which can cause severe rotor whirl in 
axial flow compressors and turbines. One is due to the circumferential variation 
of static pressure action on the cylindrical surface of a rotor, particularly with- 
in labyrinth seals. The other is due to circumferential variation of blade-tip 
clearance (the topic of this paper) . Alford presented both seal deflection criteria 
and bladed-disk torque deflection criteria for use as design guides for a stable 
rotor system. 

Investigations of rotordynamic problems in steam turbines have shown that un- 
balanced moments on wheels and shafts, caused by steam forces, are a source for non- 
synchronous rotor whirl. Steam whirl was first described by Thomas [6] in 1958 and 
is similar to the phenomenon hypothesized by Alford. Winter [?] performed model 
tests in 1968. He used ball bearings in his test rig to eliminate the shaft whip- 
ping due to hydrodynamic bearings, and applied external damping elements to the 
bearings since ball bearings have little Inherent damping. Winter proved the pre- 
sence of a load-dependent instability but was not able to quantify his results. 


THE TEST APPARATUS 


The test apparatus used in the experiment consists of a flexibly mounted b laded 

rotor driven by a variable speed electric motor, and a movable shroud (see Eigure 

1). The blade 0. D. is 6”. The motor is rated 1/7 HP at 10,000 rpm. 

An offset of the bladed rotor within the shroud clearance circle is produced by 

translating the shroud laterally. The shroud is mounted on slides, and is trans- 
lated laterally by turning a threaded rod. 

The lateral stiffness of the rotor bearing support is made low enough so that 
very low levels of Alford *s force (3-10 grams) produce a measurable rotor deflec- 
tion. The rotor deflection is measured with an eddy-current proximity probe and 


261 



calibrated to lateral force in grams/mil. One of the most challenaine and tln,^ 
consumng aspects of the project was the measurement of rotor denfctforCD^^^al) 

vJh^aUo'r'nrbal/hLlr^''^ runo;" sylrolous' 

probe signal through a „e^J e^fSl^-low-pLrfiu'^r"’' accompUshed hy passing the 
of 0 J^'^TOOrrom*' T'"! l”<*apa"<lan':ly from the torque, over a range 


Rotor speed is measured by feeding a once-per-revolution 
an electronic counter , 


keyphaser signal into 


torsional strain gages calibrated to torque. ’ 


EXPERIMENTAL PROCEDURE 


For each set of data taken, a base measurement was made with the shroud -tn rV.o 
irtWs^poiuion^^the rjtor*"^^ lateral eccentricity nominally zero). With the shroud 
ments were made ^f rotor deflection. “speed' 

roto^SS€“-°^ 

eccentricity in the negallvr(oplosile) dlrectloL''^^ repeated for an equal shroud 

strvI“lv“T-l^^ 


RESULTS 


The results of the experiments are shown in Figures 2 throm>h 1 rh^ i -i 
connecting the points in the figures were derived uLng a ll^qua^s m . 

itv ^Thrnn/ 1!°'' measured Alford force varied with rotor eccentric- 

eccentricity ias\mpl«ld TilnreTt "'d direction that the 

S L'l4u e^L thrSi^cu::: the dlffetL'rhS:Ln 

g es the direction that the eccentricity was Imparted. The efficiency 
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factor g was calculated using Alford’s equation which gives the aerodynamic force as 
proportional to a radial deflection which is normal to the force. The constant of 
proportionality is called cross-coupled stiffness. Thus we have 


F = = K 

y DpH ^ 


yx‘ 


X 


( 1 ) 


The cross-coupled stiffness is seen to be Ky^ = 


li 

DpH- 


Equation (1) implies that 3 remains constant for varying eccentricity. Figure 
5 confirms this. Therefore it was felt that the point on Figure 4 corresponding to 
20 mils eccentricity was not representative, and should not be included in the least 
squares fit. 


Figures 6 and 7 show how Alford’s force and 3 vary with torque. Again 3 was 
calculated using equation (1), Figure 8 shows how the measured Alford force varied 
with speed holding the rotor torque constant. Figure 9 was derived from figures 6 
and 8. The curve corresponding to 5000 rpm was derived using a least squares fit to 
the data from all 3 tests in Figure 6. The four remaining curves on Figure 9 were 
derived from this one curve and Figure 8, using extrapolation techniques. Figures 
10, 11, 12 and 13 show how the measured Alford force varied with rotor speed and 
torque. These figures help support the empirical equation of Alford’s force to be 
presented below. They show how the developed rotor torque varies with the speed. 


AN IMPROVED PREDICTION MODEL 


One of the remarkable features of Alford’s theory is its simplicity. Figures 
2, 3, and 6 verify the predicted linearity of Alford’s force with rotor eccentricity 
and stage torque. However, Figure 6 shows that Alford’s force does not appear at all 
^ certain level of stage torque is reached: then it suddenly appears and in- 

creases steeply with torque, Alford mentions, in reference [l], that ’’large power 
inputs to the compressor rotor appear to increase the hazard of whirl. The vibration 
problem was encountered only at the full 100 percent power rating of the engine.” 
Figure 6 is compatible with this observation. 


Figures 7 and 8 also reveal certain limitations of equation (1) , principally 
that the efficiency factor 3 actually varies with the aerodynamic load torque and 
that the developed Alford’s force is speed— dependent as well as torque— dependent. 
Also, extrapolation of the curves for Alford’s force versus torque (e.g. Figure 6) 
shows that the y-intercept is not zero as equation (1) predicts. 


To obtain an Improved prediction model, the speed must be included as a perti- 
nent variable. Dimensional analysis shows that if the speed is pertinent, then the 
inlet velocity is also pertinent. In dimensionless groups, the prediction equation 
takes the form 


FyDp 


= f 


T 



( 2 ) 
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where N = speed, rpm 

V = inlet velocity, fps 

and the other variables are as previously defined. 


j experiments reported here, measurements of the inlet velocity V werp nnr 


Equation (1) is modified to include a y-intercept, which gives 


F 


y 


xgx 

DpH 


+ C 


(3) 


where the constant C is evaluated as C = -17.37 from Figure 9. 

In addition. Figure 9 shows that the slope of the curve is affected by the rotor 


3x 

^ = -10.54 X 10"3 (N) + 89.33 


(4) 


Substituting equation (4) into equation (3) gives 
Fy = (-10.54 X 10“3 (N) + 89.33) t - 17.37 


(5) 


^"-*0:875 ?rglveJ°™' Pp - la., and 

Py = )(-1.12(N) + 9526.42) t - 1852.35 » ^ (6) 

} DpH 

where 

Fy “ Alford force, grams 
T = developed rotor torque, in-lb 
N = rotor speed, rpm 
X = eccentricity, in 
Dp = blade diameter, in 
H = bucket height, in 

If the units of Alford force are taken to be pounds (lb), the equation becomes 
Fy = {(-2.47 X 10-3(rpm-l) x (N) + 21.0) x - 4. 08 (in-lb) | ( 7 ) 

rhar ^52 ^^easonably approximates all of the curves presented in this paper 

aeca„?^Ultf ”t -tor spaJd.'L 

the deeortJed test rig. '' '’““p" ““ith fits the measured data for 


264 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE RESEARCH 


Due to the size of the test rig and the lack of measurement capability for aero- 
dynamic velocities and pressures, the conclusions which can be made from this re- 
search are quite limited in scope* They are as follows: 

1. The aerodynamic force postulated by Alford for axial flow turbomachinery 
does exist, and can be measured experimentally. 

2* Alford's force is speed-dependent as well as torque-dependent, and is 
probably also a function of the gas velocity entering the stage. 

3. Although the relationship of Alford's force to the stage torque is fairly 
linear at a given speed, the y— intercept is not zero. That is, no aerodynamic force 
is generated until a certain level of stage torque is reached. 

The authors believe that these conclusions are significant, since rotordynamics 
engineers have used Alford's theory for years with no experimental verification. 

However, it is easy to see why rotordynamic stability thresholds cannot be 
reliably predicted at present by computer analyses which use Alford's equation as a 
model for the destabilizing coefficients. In order to do this, it will be necessary 
to develop more accurate prediction models for the aerodynamic force, which will re- 
quire experimental measurements on full scale turbomachinery components, such as an 
axial flow stage from the compressor in a turbojet aircraft engine, or an industrial 
process compressor. 
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Figure 10. - Alford force as functi 
speed, showing ±2 standard deviat 



Figure IK - Alford force as functio 
speed, showing ±2 standard deviatl 
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